Rotations in isospace: a doorway to the understanding 
of neutron-proton superfluidity in N=Z nuclei 
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The T=2 excitations in even-even N=Z nuclei are calculated within the isospin cranked mean-field 
approach. The response of pairing correlations to rotation in isospace is investigated. It is shown 
that whereas the isovector pairing rather modestly modifies the single-particle moment of inertia in 
isospace, the isoscalar pairing strongly reduces its value. This reduction of the moments of inertia 
in isospace with respect to its rigid body value is a strong indicator of collective isoscalar pairing 
correlations. Beautiful analogies between the role of isovector pairing for the case of spatial rotations 
and the role of isoscalar pairing for the case of iso-rotations are underlined. 
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The ground state of most nuclei can be characterized 
in terms of a superfluid condensate of Cooper pairs which 
are formed by nucleons moving in time reversed orbits [Q . 
Nuclei with identical number of protons and neutrons, 
N=Z, exhibit an additional symmetry, related to the 
similarity of proton and neutron wave functions at the 
Fermi surface. In such a case protons and neutrons occu- 
pying identical spatial orbitals may form two fundamen- 
tally different Cooper pairs of either isovector (T=l) or 
isoscalar (T=0) type 0. The question whether isoscalar 
pairing may form a condensate similar to the well estab- 
lished isovector pairing has gained considerable interest 
in recent time. 

Already early on it was noticed that the understanding 
of excitation energies of the isobaric analogue states yield 
important information on the effective nuclear force, in 
particular also on its pairing component |3j^|] , see also re- 
cent works HD- Therefore, we analyze the T~2 excita- 
tions in even-even N=Z nuclei by means of the cranking 
approximation in isospace which is the lowest (linear) or- 
der approximation to the projection onto good isospin |Q . 
This approximation has been tested within an exactly 
solvable model by Chen et al. M where it was concluded 
that, in conjunction with number-projection, it offers a 
reliable approximation to the exact solutions. Another 
motivation to apply the cranking approximation stems 
from the formal analogy between spatial and isospin ro- 
tations following arl(l + 1) and cltT(T + 1) patterns, 
respectively ||. The crucial quantity of our investiga- 
tion is the inertia parameter in isospace, ax [reciprocal 
of the moment of inertia Indeed, the study of the 

rotational spectra was crucial to establish evidence for 
superfluidity in atomic nuclei JlcQ - Similarly, our micro- 
scopic calculations of the nuclear inertia parameter or in 
isospace show that it strongly depends on the short-range 
pairing correlations. However, in contrast to spatial rota- 
tions it is shown that ax is rather insensitive to isovector 
but extremely sensitive to isoscalar proton-neutron pair- 



ing correlations. 
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FIG. 1. The single-particle routhians (upper panel) versus 
the iso-cranking frequency for the equidistant level model. 
Solid (dashed) lines depict isovector (isoscalar) sp states, re- 
spectively. At each crossing frequency (indicated by arrows) 
the configuration changes, and hence excitation energy and 
iso-alignment (lower panel). 

Before entering the details of our model, let us con- 
sider a single-particle (sp) routhian H u = H sp — u>t x . For 
simplicity, let us also assume that the spectrum of H sp 
is equidistant [e, — iSe] and iso-symmetric. Hence, at 
oj = each eigenstate of H u is four-fold degenerate. The 
cranking term, —ut x , lifts the isospin but not Kramers 



degeneracy resulting in isovector \iv) 



s/2 



(\n)-\p)) and 



isoscalar \is) = -j^(\ n ) + \p)) doublets, with isospin align- 
ment of Tl/2, respectively (Fig. 1). Clearly, the ground 
state configuration changes stepwise at the crossing fre- 
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quencies: = 8e, 3£e, 5Se, • • • , (2n — l)<$e. At each 

crossing frequency two isoscalar states are emptied and 
two isovector become occupied. Hence, the total isospin 
alignment, (i x ) = T x , changes in steps of AT X = 2. Since 
at the same time T z = T y = 0, i.e. AT X = AT, the 
stepwise increase corresponds to non-collective rotation. 
The ground state band (gsb) consists of only even isospin 
states, T = 0, 2, 4, • • • , 2n in complete analogy to the even 
spin sequence in the gsb of spatially rotating even-even 
nuclei. Exploring further this analogy, one can show that 
the even-T sequence in the gsb is a consequence of iso- 
signature symmetry, R T = exp(— iirt x ), similar to signa- 
ture symmetry conservation, R = cxp(— iirj x ), for rota- 
tional motion. It is of importance to underline that odd- 
T states can be reached only by the proper particle-hole 
excitation at lo = 0. 

Once the crossing frequencies are calculated, it is 
straightforward to compute the excitation energy Et 
(with respect to the gs) spectrum of the iso-rotational 
gsb band: 



E T = E UJ + cuT x 
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The schematic sp model leads to the classical rotational 
formula with an inertia parameter proportional to 
the single-particle splitting at the Fermi energy. 
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FIG. 2. Alignment versus frequency (a) and excitation 
energy versus alignment (b) for the sp model (discrete 
step- function and *) and for a model including isovector pair- 
ing correlations. Calculations have been performed for 48 Cr 
at fixed deformation = 0.25. 



Let us now briefly investigate how the sp model is af- 
fected by the presence of isovector pairing correlations. 
To study this issue we have performed a series of Lipkin- 
Nogami calculations for selected N=Z nuclei using the 
deformed Woods-Saxon potential as a mean-field (at 
fixed deformation) and standard isovector seniority-type 
pairing interaction [ pd| . A representative example reflect- 
ing the generic results of the study is illustrated in Fig. 2. 
The major modification introduced by isovector pairing 
correlations is the smooth increase of iso-alignment with 
cranking frequency, see Fig. 2a. The isovector pairing in- 
troduces a kind of collectivity on the top of the sp model 
but does not affect the bulk properties. The value of the 
inertia parameter depends on the shell-structure at the 
Fermi energy but in a modest and intuitively understand- 
able way. Namely, as compared to the sp estimate, it de- 
creases (increases) when shell gaps are present (absent). 
It is also interesting to note that the smoothing effect de- 
pend only weakly on the strength of isovector pairing cor- 
relations, Gt=i- The alignment, as well as the excitation 
energy almost do not change, even when the isovector 
pairing strength, Gt=i, is strongly reduced. Henceforth, 
the excitation energy of the T-states will be computed at 
the standard cranking constraint T x = \/T(T + 1). 

The main objective of our study is to clarify the role 
played by isoscalar pairing correlations, in particular we 
show that one simultaneously: (i) can recover the Wigner 
energy as it was shown in our previous work ]Tl]J[^ ] and 
(ii) determine the excitation energies of the T=2 states 
in even-even N=Z nuclei. In addition, our subsequent 
publication will show that the mean-field model incor- 
porating isoscalar pairing correlations also is able to ex- 
plain on the same footing the excitation energies of T=l 
states in even-even N=Z nuclei as well as the compe- 
tition of T=0 and T=l states in odd-odd N=Z nuclei. 
Our Hamiltonian is based on the deformed mean-field 
potential of Woods-Saxon (WS) type JT|. The two body 
residual interaction contains both isovector and isoscalar 
seniority pairing: 

H u = hws + Gt=i P}Pi + G T =oP^Po - ut x (2) 

where P\ and Pq create isovector and isoscalar pairs, 
respectively. The Hamiltonian (0) is solved using the 
Lipkin-Nogami method. The model is very similar to 
the one described in detail in Ref. [O. However, differ- 
ent to Ref. |L2| we now employ the most general Bogoli- 
ubov transformation. It allows us to fully explore the 
isoscalar pairing channel without any symmetry induced 
restrictions i.e. to include simultaneously olol and act 
isoscalar pairs. In the present study, where we confine 
to 1=0 states in even-even nuclei, it is sufficient to con- 
sider aa T=0 correlations. Moreover, since this study 
aims at a qualitative description, we have assumed near 
spherical deformation, 02 = 0.05, for all nuclei. A draw- 
back of our model is the lacking response of the isovec- 
tor particle-hole (p-h) field to rotations in isospace. To 
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fully investigate in a quantitative way the interplay of 
isoscalar pairing and rotations in isospace, requires self- 
consistent Hartree-Fock-Bogoliubov calculations with re- 
alistic isovector and isoscalar two-body interactions in 
both p-h and pairing channels. This, however, is clearly 
beyond our present approach. 
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FIG. 3. Calculated Wigner energy strength VK(A) (a) in 
N—Z nuclei. The insert shows the isovector and isoscalar 
pairing strength parameters used to reproduce the smooth 
empirical trend of the Wigner energy strength 47 /A MeV. The 
lower panel (b) shows the excitation energy of the T=2, 7=0 
states. The experimental data values are marked by (*); the 
calculations including both T=0 and T=l pairing are marked 
with (•) and the calculations with T=l pairing only are la- 
belled with (o). 



shown in the insert. The calculated and empirical excita- 
tion energies of the lowest T—2 states, AEt=2, for even- 
even N=Z, 20<A<56, nuclei are displayed in Fig. 3b. 
Calculations including isoscalar pairing correlations (•) 
are in excellent agreement with the empirical data (*). 
In contrast, calculations including only isovector pairing 
field (o) account for roughly half of the empirical excita- 
tion energy. It is interesting to notice that although gen- 
erally AEt=2 decreases as a function of A, it clearly rises 
at closed (sub-)shells, particularly for N=Z=20. This 
nicely reflects the dominant role played by the sp sub- 
structure AEt=2 oc Se [see Eq. ([!])] and the reduced role 
of the pairing correlations (in particular the T=0 pair- 
ing) at shell closure. Note also the pronounced smooth- 
ing effect of isoscalar pairing on the AEx=2 excitations at 
shell-closure. Last note least, it is important to remem- 
ber that the Wigner energy and the T=2 excitations are 
calculated in a totally different manner. 
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The isovector pairing strength, Gt=i, is computed 
using the average gap method of Ref. Km where the 
number of proton (and neutron) WS states retained for 
the pairing calculations is consistently put to A/2. To 
compute the strength of the isoscalar pairing correla- 
tions, Gt=o, we follow the prescription given in Ref. p2[ . 
This method is based on the assumption that, within 
the mean-field model, the Wigner energy is predomi- 
nantly due to the T—0 pairing correlations. In other 
words we fit Gt=o to reproduce roughly the Wigner 
energy strength W (A)m4:7 / AMeV using the technique 
provided in Ref. |15|. The result of these calcula- 
tions is shown in Fig. 3a. Interestingly, to obtain the 
proper value of W(j4), the mass scaling of isovector and 
isoscalar pairing strengths has to be different. The ratio 
x t=0 =Gt=o/Gt=i necessary to reproduce the empirical 
trends decreases smoothly from x T=0 «1.65 at the begin- 
ning of the sd-shcll to x T=o txlA0 in the f 7 / 2 subshell as 



FIG. 4. Alignment (a) and isoscalar and isovector gap 
parameters (b) versus iso-cranking frequency calculated for 
24 Mg (•) and 48 Cr (o). The fi gure illustrates the phase tran- 
sition leading to the disappearance of isoscalar T—0 pairing 
correlations (lower panel) once the alignment reaches a value 
of \/6 (upper panel). 

There is an appealing correspondence between the role 
of isoscalar and isovector correlations. The binding en- 
ergy of even-even nuclei is lowered with respect to their 
odd neighbours due to the well known blocking effect. 
Analogous, the generalized blocking of isoscalar pairing 
correlations lowers the ground state of even-even N=Z 
nuclei and accounts for the missing binding energy com- 
monly known as the Wigner energy JTT[ . Isovector pair- 
ing is weakened at high angular momenta due to the 
Coriolis effect, that tends to align the angular momenta 
of the nucleons along the rotational axis. Similarly, be- 



3 



cause the T=0 pairs have isospins coupled antiparallel, 
rotations in isospace tend to destroy these correlations. 
In contrast, isovector pairs have their isospins coupled 
parallel and are hardly affected by iso-rotations. Pair- 
ing correlations as a function of rotational frequency in 
either space or isospace are quenched in a similar fash- 
ion like the magnetic-field destroys the electronic Cooper 
pairs in metallic superconductors. Hence, with increas- 
ing iso-cranking frequency one does expect a bulk phase 
transition similar to the well known Meissner effect [jl6| . 

The phase transition discussed above on general 
grounds indeed occurs systematically in our calculations. 
Two representative cases are depicted in Fig. ||. The 
phase transition (cf. lower panel of Fig. 4) takes place 
almost exactly at, or just before the iso-alignment (cf. 
upper panel of Fig. 4) reaches the value of T x = y/6, cor- 
responding to the T—2 state. Evidently, once we reach 
the T=2 states, isoscalar pairing correlations have es- 
sentially dropped to zero. At the same time, isovector 
pairing correlations are still strong. The quenching of 
isoscalar pairing in the ground state of \N-Z\=4 nuclei 
is a general feature of most known calculations, inde- 
pendent on the interaction |0,E^,E^ . These findings are 
therefore consistent with the isobaric symmetry which 
demands that the structure (and hence, also the exci- 
tation energy) of T=2,T Z =0 states to be similar to the 
structure of T=2,T Z =±2 members of the T—2 quintu- 
plet. The T=2,T Z =±2 states are just the ground states of 
the \N-Z\=4 nuclei, since, by the rule, the ground states 
of even-even nuclei are the states of minimum isospin: 
T=\T z \=\N-Z\/2. 



and N-Z for a sequence of Cr-isotopes is further illus- 
trated in Fig. ||. At small frequencies the T=0 pairing 
phase is present only in the T Z =0,1 isotopes, where it 
strongly lowers the Mol. At high frequency a phase tran- 
sition takes place resulting in a rapid increase of the Mol 
to the 'rigid body' value i.e. the value corresponding to 
S*jf ^ for the T=0 unpaired system. 

In summary, the response of pairing correlations to ro- 
tations in isospace is investigated within a simple model. 
The present calculations show that on a qualitative level, 
the mean-field method is capable to account for both 
mass-excess in N—Z nuclei and the Mol in isospace if 
and only if the short range correlations take into account 
isoscalar pairing. Pairing correlations of isovector and 
isoscalar type respond totally different to rotations in 
isospace. The presence of isoscalar pairing strongly re- 
duces the Mol in isospace, but only for low values of T. 
With increasing iso-cranking frequency, isospin starts to 
align, iso-pairs become broken, resulting eventually in the 
quenching of isoscalar pairing and a rapid increase of the 
Mol. In the regime of large isospin, no isoscalar pairing is 
present. These results are in beautiful analogy to spatial 
rotations. 

This work was supported by the Goran Gustafsson 
Foundation, the Swedish Natural Science Research coun- 
cil (NFR), and the Polish Committee for Scientific Re- 
search (KBN) under Contract No. 2 P03B 040 14. 
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FIG. 5. Moment of inertia (Mol), 3;, ee T x /lu, for a se- 
quence of N-Z— 0,2,4,6,8 Cr-isotopes at constant deformation 
(/?2=0.25) versus rotational frequency. The figure clearly il- 
lustrates the pronounced effect of the T=0 pair correlations 
in N-Z—0 and 2 isotopes and the phase transition rising the 
Mol to the 'rigid-body' value at high frequencies. 
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The dependence of the moments of inertia (Mol), 
T x /lo as a function of the iso-cranking frequency 
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